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Abstract 

We consider spin accumulation at a ferromagnet-normal metal interface in the presence 
of magnetic scattering in the normal metal. In the classical regime, we discuss the inverse 
Drude scaling of the conductance as a function of the interface transparencies. We present 
a treatment based on an exact solution of the Boltzmann equation. In the quantum 
regime, we solve a single impurity "spin-flip Fabry Perot interferometer" for quantum 
coherent multiple scatterings, in which we find a resonance in the spin flip channels. This 
resonance appears to be the quantum analog of the semi classical inverse Drude scaling 
of the conductance. 
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1 Introduction 



The discovery of the Giant Magneto Resistance (GMR) in magnetic multilayers |T|, ^, ^ has generated 
an important interest. These systems are made of a sandwich of alternating ferromagnetic and non 
magnetic layers. Valet and Fert proposed a semi classical description of the perpendicular GMR, 
on the basis of a Boltzmann equation incorporating a spin-dependent transport in the presence of 
spin accumulation (see also [^). Spin accumulation occurs in the GMR because the current 
arising from a ferromagnet is spin polarized, and therefore cannot penetrate a ferromagnet with an 
opposite magnetization. Instead, spin accumulates at the interface. This phenomenon occurs also 
at the interface between a ferromagnet and a superconductor, where a spin polarized current cannot 
penetrate the superconductor |6|, |^, ^. Here, we would like to reconsider two particular aspects of spin 
accumulation, namely, (i) in the semi classical regime, the possibility of an inverse Drude scaling of 
the conductance meaning that, in some parameter range, the conductance increases with the length 
of the conductor; and (ii) in the quantum coherent regime, the existence of a resonance in the spin 
flip channels. More precisely, we study a ferromagnet - normal metal - ferromagnet spin valve, in 
which we assume the presence of magnetic scattering in the normal metal j^, |l^] . The inverse Drude 
scaling resulting from spin accumulation is already implicitly contained in the equations obtained by 
Valet and Fert but, to our knowledge, this effect has not been studied previously per se in the 
literature, which we do here. Spin accumulation corresponds to the presence of a different chemical 



potential for the spin-up and spin-down electrons, which obey a spin diffusion equation Q, 11, 12 1. 
Our treatment is not based on the spin diffusion equation, but relies on an exact solution of the ID 
Boltzmann equation where we can make an exact decoupling between the charge and spin sectors. 

Next, we ask to what extend a quantum model can show a similar physics. We are lead to study 
a quantum "spin-flip Fabry Perot" interferometer in which a single magnetic impurity is located at a 
given distance a away from a ferromagnet interface. We find the existence of a Fabry-Perot resonance 
in the spin-flip channels as the parameter a is varied. This resonance disappears as the ferromagnet 
spin polarization is decreased, and can therefore be viewed as the equivalent of the inverse Drude 
behavior in the quantum coherent regime. Our treatment is based on a Landauer approach, similar 
to the one used by Zhu and Wang to study the effect of magnetic scattering close to a superconductor 



interface [|13|. 



The article is organized as follows. Section ^ is devoted to the solution of the semi classical 
transport equations. We solve the "spin-flip Fabry Perot" interferometer model in section ^. Final 
remarks are given in the Conclusion. 

2 Transport in the semi classical regime 
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Figure 1: The "spin valve" geometry considered in section |2|, consisting of a normal metal wire doped 
with magnetic impurities, connected to two ferromagnets with an opposite magnetization. 



2.1 Boltzmann equation and boundary conditions 

We consider a model in which magnetic impurities are present in a normal metal close to a ferromagnet 
interface (see Fig. ||). We neglect any Kondo correlation |14], which is an assumption valid above 
the Kondo temperature. The presence of the ferromagnets close to the normal metal may lead to 
magnetic flux lines penetrating inside the normal metal, which can orient the magnetic impurities 
in a preferential direction. We implicitly assume that the temperature is high enough so that the 
impurities have no preferential orientation. We consider a one dimensional model because only in 
this geometry can we decouple the spin and charge sectors of the Boltzmann equation. We assume 
that the interfaces between the ferromagnets and the normal metal are sharp, and that the exchange 
field has a step function variation at the interface. We note f^j^{E,x) the semi classical distribution 
function of right/left moving spin-o" electrons with an energy E at position x. 
The Boltzmann equation in the relaxation time approximation reads 
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where we have discarded the term involving the electric field. This is valid if the temperature of 
the electrodes is larger than the applied voltage, in which case the electronic gas has a temperature 
identical to the one of the electrodes 16 1. Therefore, we should consider a finite temperature 



and calculate the low voltage conductance in the regime eV <^ T. In practise, we consider the limit 
T — > 0, and calculate the linear conductance. The coefficients r, and r'^ in Eq. [l| denote respectively 
the rate of backscattering without spin- flip, the rate of forward scattering with spin- flip, and the rate 
of backward scattering with spin-flip. The coefficients can be related to the q = and q = 2kf 
components of the microscopic scattering potential (see the Appendix). 
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We now explicit the boundary conditions. For this purpose, let us consider an interface between 
a ferromagnet in the region x < and a normal metal in the region x > 0, and include interface 
scattering under the form of repulsive potential H6{x) Let us first consider a spin-up electron 
incoming from the left ferromagnet, and denote by and the backscattering and transmission 
coefficients. The wave function in the region x < is ipLix) = exp (ik^x) + exp {—ik^x), and the 
wave function in the region x > is iPr{x) = t^exp{ikx). The matching equations are V'-l(O) = 
i^niO) = tp{0) and dil)R{Q)/dx — dil)L{^)/dx = (2mff/7i^)^(0), from what we deduce 

V = ^, and W = 

i{k + k^) -2mH/n^ i{k + k^) -2mH/f? 

The probability current conservation can be verified easily: k^ = k^\b^\'^ + k\t^\'^. The spin-up 
conductance is found to be = {e^ / h)\T'^ , with the transmission coefficient 



(A; + A;T)2 + 



2mH 



The backscattering coefficient is i?^ = 1 — T^. In the spin-down sector, we obtain and by 
substituting k^ with k^ in Eq. ^. This provides the boundary conditions for the Boltzmann equation: 

fl[E,Q) = T^fj.(E-eV) + il-T^)fliE,0) (3) 

fj,{E,0) = TUT{E-eV) + {l-T^)fi{E,Q) (4) 

fl{E,L) = T^fT{E) + {l-T^)fl{E,L) (5) 

fi{E,L) = T^fT(,E) + {l-T^)f^^{E,L), (6) 

where the left and right ferromagnets are assumed to be in equilibrium and fxiE) denotes the Fermi- 
Dirac distribution function. We consider T^^ and to be independent of energy, which amounts to 
considering the wave vectors k and k^ in Eq. ^ to be on the Fermi surface. Eqs. ^ ^ ^ provide a 
simple form for the spin-up and spin-down currents at positions x = 0,L. For instance at x = L, we 
have 

iHl) = T^^J[fl{E,L)-fT{E)]dE (7) 

I^{L) = T^^J[fj,{E,L)-fT{E)]dE. (8) 

The Boltzmann equation Eq. |l| and the boundary conditions Eqs. ^ |^ lead to eight equations for 
eight variables fj^^^{E,x = 0,L). We now solve these equations directly and discuss their physics. 

2.2 Solution of the Boltzmann equation 

The 4x4 Boltzmann equation can block diagonalized into 2x2 blocks by changing variables to 

T I T I 

the charge and spin combinations Xr^l = fRL~^fRL Yr,l = Ir l ~ fk l- This spin-charge 
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decoupling allows to solve exactly the Boltzmann equation. We find 
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with I = l/{r + r'g) the mean free path, and a = r + 2r<j + r^, b = r — r'^. The 2x2 block equations 
can be easily integrated to obtain 
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with X = L/l. Similarly, 
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where t = cosh (AL) — asinh(AL), t = cosh (AL) + asinh(AL), and u = /3sinh(AL). We used the 
notation a = a/ \, l3 = b/\, and A = Va^ — 6^ . Next, we combine the boundary conditions Eqs. 
to Eq. |10| to obtain an expression for fj^{E,0) — fj^{E,0) and fi{E,0) — fl{E,0) as a function of 
fj^{E,L) and f^{E,L). Once injected into Eq. these relations lead to 



2ME-eV)f( T^ + Tl-2T^Tl 



(12) 



The matrix M appearing in the left hand side of Eq. ^ is 



M = T 



(T^ - T^) 



-1 1 



with the coefficients 
A^ 



3TT _ 4tTt^ + - 2(tT)2 + 2(TT)2ri + xT^ [t^ - 2T^T^ + 



(13) 



(14) 
(15) 



The expression of A^ is obtained by exchanging and in Eq. 14. Similarly, B^ is obtained from 
B^ by exchanging T'^ and in Eq. |l^. 

2.3 Fully polarized limit 



We first consider the solution Eq. |12| in the case of fully polarized ferromagnets with high transparency 
contacts: H = = = k} = k, leading to = 1. In this limit, only a spin-up current 



4 



r=0.1;rs=0.01 



(U 

■t— ' 

"o 

> 

a 
'■*—> 

N 



o 

Oh 



0.25 




Figure 2: Spin polarization profile inside the diffusive conductor, with r = 0.1, = 0.01, and 
increasing values of L. The spin polarization is normalized to the applied voltage, and the coordinate 
along the wire is normalized to the total length. We have 1/A = 15.8. 



can enter the ferromagnet at x = L. This is expected on physical grounds, and it can be verified 
explicitly on the form Eq. ^ of the spin-down current. The total current is found to be 
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(16) 



where we considered only the forward scattering spin flip processes 



0), and assumed that 



<C r, in which case the elastic mean free path ^ = 1/r is much below the spin- flip length /gf 



l/[2-\/rs(rs + r)]. If L is small compared to Zgf, the conductance G ~ 2^rsL shows an inverse Drude 
behavior. 



2.4 Spin polarization profile 

The spin polarization profile in the diffusive wire can be calculated in a straightforward fashion from 
the solution of the Boltzmann equation. Once we know the distribution functions at one extremity 
of the wire, we can use Eqs. ^ to propagate the solution to an arbitrary point. The resulting spin 
polarization inside the wire is proportional to the applied voltage, and is shown on Fig. ^ for various 
values of L. When L > /gf, there is a plateau in the spin polarization in the middle of the wire. In 
the opposite inverse Drude regime, there is no such plateau. 
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Figure 3: Variations of the conductance as a function of the length L of the diffusive wire, with 
decreasing spin polarizations = 0, 0.1, 0.2, 0.3, 0.4, 0.5, and the parameters = 0.9, r = 0.1 
and Vg = 0.01. With a strong spin polarization, the conductance increases with L below Lc- With a 
weak spin polarization, the conductance decreases monotonically with L (Lc = 0). 
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Figure 4: Variations of Lc as a function of T^, with increasing spin polarizations = 0.1, 0.2, 0.4, 
0.6, 0.8, 1, and the parameters r = 0.1 and = 0.01. With a given T'^ , the curve Lc{T^) separates 
two regions: (i) small-L, small-T-l-: the conductance increases with L; and (ii) large-L, large-T-l-: the 
conductance decreases with L. 
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2.5 Effect of a partial spin polarization 

Now, we consider the effect of a partial spin polarization in the ferromagnet. It is expected on physical 
grounds that a decreasing spin polarization tends to suppress the inverse Drude scaling because this 
regime is clearly absent in the spin unpolarized case. This is visible on Fig. ^ where we plotted the 
conductance as a function of the length of the diffusive wire for decreasing spin polarizations. With 
an arbitrary polarization, there exists a critical length scale Lc such that the conductance increases 
with L below Lc, and decreases with L above Lc- There exists also a critical value of such that 
Lc = if > T^. To illustrate this, we have shown on Fig. Q the variations of the critical length Lc 
as a function of the parameter TK When increases, Lc decreases: the maximum in G{L) occurs 
for a smaller Lg. When is above a critical value TJ^, the conductance decreases monotonically 
with L. 

We now describe the effect of a partial spin polarization on the basis of a small-TJ- expansion. 
The strategy is to express the current to order L and determine whether the conductance increases or 
decreases with L. We expand the current to first order in the two parameters K = XL and x = L/l, 
and retain the coefficients of this expansion to leading order in T^. It is first instructive to carry 
out the expansion with L = 0, and therefore K = x = ^. It is visible on Eqs. ^, |8| that a prefactor 

enters the spin- up current, and a prefactor enters the spin-down current. We should then 
express fji{E, L) to first order in while ffi{E, L) should be expressed to order {T^)^ . The spin-up 
and spin-down channels appear to play an asymmetric role. Nevertheless, the final expression of 
the conductance is identical in the spin-up and spin-down channels. An intermediate step in the 
calculation of fji{E, L) is the derivation of DetM to order (see Eq. 13): 



BetM = -AiT^f Y - : 

leading to an identical current in both spin channels: 

j1 = ji = I [fT[E - eV) - ME)] dE. 

Now we consider a diffusive wire with a finite length L, and expand the current to first order in x 
and K, and to leading order in . The determinant of the matrix M in Eq. 12 is found to be 

DetM = -A{T^f \l-T^ | - A{T'^f{a - (3){2 - T^)K - Ax{T^ fTK 

Next we expand the spin-up current to order L to obtain 

If is small, the current increases with L while it decreases with L if is large. The transition 
between these two behaviors is obtained for Tj- = \J (2rs/r)T^ , compatible with the behavior shown 
on Fig. H. 
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2.6 Replacement of one of the ferromagnets by a normal metal 

We now consider the situation where we replace the left-hand-side ferromagnet on Fig. |^ by a normal 
metal. In the presence of high transparency contacts, the conductance of this junction is of order 
/h in the absence of diffusion while it is of order {e'^/h)T^ in the spin valve geometry on Fig. |l|. 
Replacing one of the ferromagnets by a normal metal is expected to suppress the inverse Drude 
scaling. The boundary conditions appropriate to describe this situation are 

/^(^,0) = TfT{E-eV) + {l-T)fl{E,0) (17) 

/^(S,0) = TfT{E-eV) + {l-T)fl{E,0) (18) 

fliE,L) = T^fT{E) + {l-T^)fj,{E,0) (19) 

fi{E,L) = T^ME) + {l-T^)fj,{E,0), (20) 

that should be solved together with Eqs. ^. The solution is found to be 

[fkiE,L))-[{T^-Ti)iu + til-T)) )fT{E) + 2T^^jME-eV), 



with 



and 




= T + T^ - TT^ + xTT"^ (21) 
D'^ = t-u{l-T)-{l-T'^){t{l-T)+u). (22) 

We have plotted on Fig. |5| the conductance of the junction with high transparency contacts, where 
it is visible that the conductance decreases monotonically with the length of the diffusive wire. 

Now, reducing the contact transparency restores a regime in which the conductance increases 
with the size of the diffusive wire. This is visible on Fig. ^ where we used T = 0.01 and = 1. 
Again, the inverse Drude scaling is obtained for the smallest values of (with strongly polarized 
magnets). 

3 Quantum coherent transport: a single magnetic impurity "spin- 
flip Fabry Perot interferometer" 

3.1 Matching equations 

We now consider a single magnetic impurity at x = in a normal metal, in the presence of a normal 
metal - ferromagnet interface at x = a (see Fig. |^. The purpose of this calculation is to study a 
model in which the interplay between multiple reflections and phase coherence is treated exactly, and 
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Figure 5: Conductance of the junction with a diffusive wire connected to a normal metal and 
a ferromagnet, with high transparency contacts: T = = 1. It is visible that the conductance 
decreases monotonically with the length of the diffusive wire. 
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Figure 6: Conductance of the junction with a diffusive wire connected to a normal metal and a 
ferromagnet, with low transparency contacts: T = 0.01, and = 1. A regime with a conductance 
increasing with the length of the diffusive wire is visible for small values of the parameter TK 
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Figure 7: The system considered in section ^. The impurity is in the normal metal at a distance 
a away from the ferromagnet. We have represented a quasi one dimensional geometry while the 
calculation is made in a one dimensional geometry. 



to determine whether there exists a signature of spin accumulation in the quantum coherent regime. 
We find that the quantum model behaves like a Fabry Perot interferometer, with a resonance in 
the spin flip channels. This can be viewed as the signature of spin accumulation in the quantum 
coherent regime. We neglect the Kondo effect because we want to describe a situation in which the 
temperature is above the Kondo temperature. The conduction electrons are scattered through the 
Hamiltonian H = Vb + ViSi.s, where Si is the impurity spin and s the spin of the conduction electron, 
and we further assume a single channel geometry. This type of model has been used by Zhu and 
Wang to investigate the effect of a magnetic impurity close to a normal metal - superconductor 
interface. The spin-up and spin-down wave functions are grouped in a two-component spinor ip{x). 
Clearly, the impurity couples the spin- up and spin-down wave functions. The matching of the wave 
function at the impurity site reads 



2m r 



V'(O-), and ^(0+) - ^(0~) = ^ Xl + fxa, V'(O), 



(23) 



with A = Vo — Vi/4 and ^ = Vi/2. The matching of the wave function at the ferromagnet boundary 
reads 

V'(a+) = i,{a-), and |^(a+) - ^(a-) = '^Hi,{a), (24) 



where we included a repulsive interface potential H5{x — a) at the normal metal - ferromagnet 



interface. Eqs. 23, 24 generate eight constraints, for a set of eight transmission coefficients. 

This calculation amounts to a resummation to all orders of a series of diagrams in which a 
conduction electron scatters onto the impurity, scatters back onto the interface, scatters again onto 
the impurity, ... (see Fig. ^ (a)). Note that the diagram with a hole in the intermediate state shown 
on Fig. ^ (b) generates another series which is not included in the calculation. If one wanted to 
describe the Kondo effect close to a ferromagnet interface, it would be crucial to incorporate the 
diagram on Fig. ^ (b), as well as inserting the interface scattering in this diagram. 
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(a) 




(b) 




Figure 8: (a) The processes included in the Landauer calculation. The wavy lines indicate the 
scattering at the interface; and (b) A process with a hole in the intermediate state, not included in 
the calculation. 

3.2 Scattering in the total spin — sectors 
3.2.1 Incoming electron with a spin-up 

We first consider a spin-up electron incoming on the interface while the impurity is supposed to have 
initially a spin down. The wave functions are 



where and denote the spin-up and spin-down Fermi wave vectors in the ferromagnet. In the 
notation of the transmission coefficients, the superscript denotes the initial and final spin orientations 
of the conduction electron while the subscript denotes the initial orientation of the impurity. The 
solution of the matching equations is straightforward, and we find the transmission coefficients 




(25) 



(26) 



(27) 
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We used the notation = ^ + iZ + = ^ - (^iZ + , ^ = exp {ika), A" = exp {ik'^a). 

The dimensionless scattering potentials in Eqs. ^ are z = mX/{h'^k) and z' = m^/ifi^k) at the 
impurity site, and Z = mH / [fi^k), Z" = mH / {fi'^k'^) at the normal metal - ferromagnet interface. 
The denominator D in Eqs. ^ is 

V = X^X^ (Af [l-z^ + {z'f + 2iz\ + (xTyJ- + X^Y^)[-z'^ + {z'f + iz] + F^y^ [-z^ + (z')^] • (30) 

3.2.2 Incoming electron with a spin-down 

We consider now an incoming electron with a spin-down while the impurity has initially a spin-up. 
The wave functions are 











"it 


:;) 



= I e--+ e--'ifx<0. (31) 



V^^l(x) = ^ je^^-+(^ ^' je-^-ifO<x<a. (32) 



The equations for t^l^*^^ and are obtained from the ones in section 3.2.1 under the transfor- 



mation ^ A^, and Z'^ Z^-. The amplitude for transmission in the ferromagnet is 

ttt""^ = -:^^z'[AX^+AY^] (34) 
t-l^'^ = :^[A{1 + iz)X^ + MzY^]. (35) 

3.3 Scattering in the total spin = ±1 sectors 

The incoming electron does not undergo spin- flip scattering in the sectors with a total spin = ±1. 
The transmission coefficients in the sector 5^ = 1 is found to be 

*T AT lA{l + i{z + z')) XT +Ai{z + z')Y^ ' 

In the sector = —1, we have 

fl^^l = I 

Al [A{l + i(z + z')) Xi +Ai(z + z')Yi] ■ 



We can check easily that these forms of the transmission coefficients are identical to Eq. 28, with 
z' = 0, and the replacement z ^ z + z'. This is expected since there is no spin-dependent scattering 
in the limit z' = of Eq. p^. 
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Figure 9: Spin-flip, non spin-flip, and total conductances (in units of e^/h) of the quantum mechan- 
ical model of magnetic scattering close to a ferromagnet interface, with an unpolarized ferromagnet 
{k = = = 1) and the parameters Z = 1, a = 100 and Wq = Vq/H = 0. The conductances are 
plotted as a function of Wi = V\/H. 

3.4 Landauer formula 

We now evaluate the total conductance and assume that the incoming electron and impurity do 
not have any preferential direction. The conductance is the sum of four terms, weighted by the 
probability V = 1/2 to have a spin-up or spin-down impurity: G = -|- + + /2, 
with 



h [ k 



.eT^eT|2 



k\ 



h 




h k^ 

e2 Refci 
¥ k 



1*: 



ej.^ei|2 



(36) 

(37) 

(38) 
(39) 



We have incorporated the possibility of having a pure imaginary wave vector fe^-, corresponding to 
an empty spin-down band. 



3.5 Resonances 

We consider the presence of a strong interface scattering at the metal - ferromagnet interface. The 
electron are multiply reflected before they enter the ferromagnet, and therefore the resonator has 
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Figure 10: Spin-flip, non spin-flip, and total conductances (in units of e^/h) of the quantum 
mechanical model of magnetic scattering close to a ferromagnet interface, with a strongly polarized 
ferromagnet {k = = 1, k^ = 0.01) and the parameters Z = 1, a = 100 and Wq = Vq/H = 0. The 
conductances are plotted as a function of Wi =Vi/H. A resonance, not present on Fig. ^, develops 
in the conductance upon spin polarizing the ferromagnet. 



a high quality factor. We first choose the parameter a in such a way that spin-flip scattering is 
resonant. As it is visible on Figs. ^ and IC, the presence of a spin polarization in the ferromagnet 
generates a resonance in the conductance, not present in the unpolarized situation. We have shown 
the conductance with Vq = but a similar behavior has been obtained with a finite Vq- The values 
of ka for which a resonance occurs can be worked out by calculating the transmission coefficients in 
the limit of a large z, z' . In this limit, we find 



1 



l\2 



izA^ 



z2) AX^ + AyT 
1 



'^2 



z2 AX^ + AY^ 



(40) 
(41) 



The resonances occur when tan (ka) = — 2Z). For a large Z, the resonances are close to the real 
axis: tan (fca) = — 1/(2Z), in agreement with Fig. |lT]. The reason why the resonance appear to be 
sharp as a function of a when Z is large is that, even without spin flip scattering, the quality factor 
of such a resonator is large when Z is large. 

The occurrence of a parameter range in which a peak occurs in the spin flip conductance is 
intriguing. The presence of a specific physics in the spin flip channels can be already understood 
from the large-z, z' behavior, Eqs. |l[ Typically, one has ~ {z')'^ /{{z')'^ - z^)^ and 



(z) /{{z') — z ) . In the presence of spin flip scattering, one has z ^ z' and therefore a 
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Figure 11: Spin flip conductance of the Fabry Perot model, as a function of a, with k = = 1, 
k^ = 0.01. As it is visible, the conductance is a periodic function of ka with a period vr. The 
resonances with Z = 1 are much sharper than with Z = 0.01, as expected on physical grounds. The 
resonances occur when tan (fca) ~ —1/{2Z). 



different conductance in the spin flip and non spin flip channels. 



4 Conclusions 

To conclude, we have determined to what extend spin accumulation can result in an inverse Drude 
behavior in a semi classical spin valve model. Our treatment was based on an exact decoupling be- 
tween the charge and spin sectors of the Boltzmann equation. We have addressed a similar question 
in a single impurity quantum model and found the existence of a resonance in the spin flip conduc- 
tance. It is an open question to determine the quantum coherent behavior of a spin valve with a 
finite concentration of impurities. 



A Derivation of the Boltzmann equation with spin-flip scattering 

We give a brief derivation of the Boltzmann equation Eq. |l| in the presence of a spin-flip scattering 
potential. The derivation generalizes Ref. |l^ to incorporate a spin-flip scattering self energy. The 
Dyson equation in the spin tensor Keldysh space reads (Gq^ - S)(l,2) ® G(2) = 6{1 - 2). The 
convolution includes a sum over coordinates, time, and spin. The kinetic equation is obtained from 
the difference of the Keldysh components of the Dyson equations and its conjugate: 



Gn^ - ReS,G 



K 



S^,ReG 



r,G 



K 



(42) 
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(a) 




(b) 




Figure 12: The self energy terms incorporated in the gradient expansion calculation. The term 
shown on Fig. ^ (b) generating the Kondo effect is not included. 

with []_ and [] , denoting a commutator and an anticommutator respectively. We refer the reader 



to Ref. pq] for an explanation of the symbols used in Eq. 42. We use the the self energy shown on 
dp' 



Fig. 12 



Scr(p,R,r) = Ui 



^;(p-p')rG^,<x(p',R,r)+n^ 



dp' 



r;'(p-p')rG_,,».(p',R,T), (43) 



(27r)3'^^ , ■^^J ^^^^ 

with rii and n'- the concentration of non magnetic and magnetic impurities. The first term in Eq. 43 
describes non spin-flip scattering, and the second term describes spin-flip scattering. Notice that this 
self energy does not incorporate the Kondo effect since we do not incorporate the possibility of a 
having hole in the intermediate state. 



We assume the self energy in Eq. ^ to be constant in space, and use the gradient expansion to 
first order 



1 + 2 [dxd. 



A R / 



If A and B are symmetric in spin A^j^^i = Aui^^^, and B^y = B^ji the commutator reduces to 
the usual spinless Poisson bracket: [A^ B]_ = i J2a2 {-^(7,02^^^2,(7}^ with {A, B} = d^AdpB — 
dpAd^B. Using these relations, we expand the kinetic equation Eq. ^and integrate over energy to 
obtain the Boltzmann equation 

dp' 



dTfp,a + Vp^pVR/p,a - VrC/Vp/p,, 



27rn,- 



p-p')\'S{^p-Cp')[U',a-fp,.] (44) 



(27r)3 
dp' 



+2nn', J ^|^;'(p-p')|25(ep-Cp,)[/p',-.-/p,<.]- 



In the one dimensional limit, Eq. ^ reduces to the Boltzmann equation Eq. |^, with the scattering 

coefficients related to the q = and q = 2k f components of the scattering potential: r = ni\v2kf\'^, 

2 



Tg = n'j\vQ\'^ , and r'g = n'^\v'2 
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